A new method for measuring the CMB temperature quadrupole with an accuracy 

better than cosmic variance 
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We present a new method for measuring the CMB temperature quadrupole, using large scale 
CMB polarization. The method exploits the fact that CMB polarization is partially sourced by the 
local temperature quadrupole. We link the temperature with the polarization spectrum directly by 
relating the local quadrupole at the onset of reionization to both of them. The dominant contribution 
is at I < 30 and since we use many I values, we can reduce the error significantly below cosmic 
variance. In particular, for our fiducial model, the error on the temperature quadrupole is reduced 
to 24%. This has the potential of reducing the probability of a low quadrupole by two orders of 
magnitude. 



How precisely can Cj, the Cosmic Microwave Back- 
ground (CMB) temperature quadrupole, be measured? 
At first sight, the question seems trivial to answer : 
the measurement is limited by cosmic variance, which 

at the quadrupole is equal to y^|Cj. The first release 
of Wilkinson Microwave Anisotropy Probe (WMAP) 
data has regenerated interest in the value of the CMB 
quadrupole, as WMAP measured a lower quadrupole 
than expected, based on a ACDM cosmology. This con- 
firmed the COsmic Background Explorer (COBE) Q 
measurement but more cleanly, as the lower detector 
noise and wider frequency range to pin down galactic 
foreground emission 0] renders the measurement more 
robust. 
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Let us consider the value of the quadrupole AXf where 
l -^p-Ci. Consider also the best fit adiabatic 
model of y as a fiducial model for the rest of the paper, a 
six parameter model with physical baryon and cold dark 
matter densities of Wf, = 0.023 and oj c = 0.117 respec- 
tively, relative cosmological constant density Cl\ — 0.715, 
optical depth to reionization r = 0.137 and scalar spec- 
tral index n = 0.974. Using frequentist statistics given 
the fiducial model and assuming full sky coverage with 
cosmic variance as the only source of error, the proba- 
bility of measuring a quadrupole as low as or lower than 
the quoted WMAP value of AXf = 123.4 {^K) 2 is 0.01. 
Using a sky cut would actually raise the probability. Ef- 
stathiou [5j has argued that statements such as the one 
above should not be taken too seriously and do not dis- 
favor the ACDM concordance model. If for example, one 
uses AT| = 201.6 (pK) 2 quoted by Tegmark et al. 
using their different method of foreground subtraction, 
the probability rises to 0.03. Efstathiou also showed that 
the value of the quadrupole is sensitive to the estimator 
used MK In particular, the use of the quadratic esti- 
mator [81 rather than the pseudo-C; estimator used by 
WMAP 9] is better suited for low I values and could fur- 
ther increase the quadrupole to ATf = 250 (fj,K) 2 . The 
corresponding probability for the fiducial model would 
then increase to 0.05, five times larger than the origi- 



nally quoted value. Other apparent coincidences such 
as the alig nment between the quadrupole and the oc- 
topole |l0( the north-south asymmetry 0] or possible 
non-Gaussianity/global anisotropy 0] do not concern us 
here even though a normal quadrupole would certainly 
weaken these findings. 

The probability arguments discussed above depend 
strongly on the cosmic variance limit. If there were a 
way to measure with precision better than cosmic 
variance, the underlying probability distribution of the 
quadrupole would be different and the probabilities could 
change drastically. Is there a different measurement one 
might make, that could measure the CMB temperature 
quadrupole and has an error smaller than cosmic vari- 
ance? This question has been raised and partially an- 
swered by Kamionkowski and Loeb . They considered 
the polarization spectrum produced by clusters through 
the Sunyaev-Zel'dovich effect 0] . Since the polarization 
spectrum depends on the local quadrupole at the clus- 
ter, one can get information on the quadrupole Cj (z) at 
redshift z, by taking large samples of clusters and taking 
an average. One hopes that the averaging method will 
recover a quadrupole close to the true cosmological value. 
The authors mention another way of getting a measure of 
, namely large angle CMB polarization, generated by 
reionization. They argue however that this would depend 
on the reionization details and will therefore be very dif- 
ficult to handle. Further calculations involving clusters 
have been carried out in and more recently in 0]. 

The connection between the temperature quadrupole 
and polarization has been exploited further by Dore et. 
al. [l8(. They considered the consistency of the observed 
temperature spectrum with the polarization-temperature 
cross correlation, given the best fit model of WMAP. In 
particular, if the observed temperature quadrupole was 
anomalously low, then one would expect low large an- 
gle polarization power as well. Their method consists of 
Monte-Carlo realizations of their fiducial model combined 
with the well known statistical correlations between the 
two above spectra [lflj . Using a frequentist approach, 
they find that the two are inconsistent at the 98.5% level. 
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This hints that the low temperature quadrupole might 
not be just due to an unlucky throw of dice. 

In this paper, we exploit the same connection between 
temperature and polarization as above and demonstrate, 
using a different method, the possibility of getting a mea- 
sure on Crf from large angle polarization generated by 
reionization. We show that contrary to the argument 
in the method does not depend on the details of 
reionization but rather on the initial amplitude of the lo- 
cal quadrupole at the onset of reionization. Any reioniza- 
tion history dependence can be accommodated by using 
extra parameters but this does not appear to be an obsta- 
cle to our method. That the details of reionization can be 
observed in the CMB has been studied before ■ Even 
in the worst case scenario, the information about reion- 
ization contained in the CMB boils down to at most five 
parameters [2l| but this is something that would have 
to be taken into account for standard CMB parameter 
extraction as well. Therefore in what follows, we use 
a sharp reionization transition at (conformal) time t ri 
the time where the electron ionization fraction first rises 
above its residual recombination value. 

We first review how the relevant spectra are related. 
Let Af(k,t) and A[(k,t) denote the /-multipole of the 
temperature and polarization transfer functions respec- 
tively. Let also t be conformal time with to being the time 
today, k the wavenumber and r(t) = J t r dt the optical 
depth to time t where a dot indicates differentiation with 
respect to conformal time. 

The temperature quadrupole is then given in terms of 
an initial power spectrum P^(k) and transfer function 
A^(fc) as 

9 r°° 

C% = - dkk 2 P^(k) |Af(fc)| 2 . (1) 
t Jo 

The quadrupole's transfer function obeys the differential 
equation 

Ar + ^rA^^(2Af-3Aj) + ^(A 3 + Ar) . (2) 

The general solution is the sum of an initial condi- 
tion (the initial quadrupole) multiplied by the relevant 
Green's function and a particular integral A| 3/ (fc) over 
the source 5*2 which is given by the RHS of © . For rea- 
sons to become clearer below, we choose the initial time 
to be the reionization time t r , so that the initial condi- 
tion is the local quadrupole at the onset of reionization. 
The local quadrupole today is then given by 

A^(Mo) = e^A^MrH^' dt e* r <*> S 2 {k,t). (3) 

Similarly the E-mode polarization spectrum is given 
in terms of the same initial power spectrum and E-mode 
transfer function Ei{k) as 

2 f°° 

Cf = - / dk k 2 P*(jfe) \Ei(k)\ 2 . (4) 
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FIG. 1: The local quadrupole today(solid) and at reioniza- 
tion(dotted). Also shown are the first five Gi's (I = 2 to 
I = 6) in dashed with increasing strength. 

The E-mode transfer function is give n by the line-of-sight 
integral along the past light cone as |22| 

where x = k(to — t), II(fc,t) is the polarization source and 
ji(x) is the spherical Bessel function. The polarization 
source is given in terms of temperature quadrupole and 
polarization monopole and quadrupole as 

n(fc, t) = A^(fc, t) + Aftfc, t) + (k, t) (6) 

and obeys the differential equation 

ri+AfII=^[2Af-3K + Af + Af)] . (7) 

The general solution of the above equation for any time 
t > t r is 

n(M) - e ral T "- T WlAj(fc,t r ) 

+e-TOH*)] J' dt ' e *r(t') s n(fcj ^ _ (g) 

where again we have chosen the initial time to be t r and 
where Su is given by the RHS of J7J|. Note that in the 
above equation we have replaced II(fc,t r ) which should 
have been the true initial condition with A^(fc,i r ). It 
turns out that due to free-streaming from recombination 
to reionization, the two are equal to one part in 10 6 . 
During reionization, however, this is no longer true as 
the three terms comprising II(fc, t) become comparable 
because of rescattering. We are therefore forced to put 
the initial condition at the onset of reionization. As one 
can see, the temperature quadrupole produced by reion- 
ization is fully connected with the polarization spectrum 
through the initial condition A^(fc,i). This forms the 
basis of our method : given a model and the polarization 
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autocorrelation spectrum Cf one can get information 
about the local quadrupole at reionization from which 
the quadrupole today can be inferred. 

One conceptual difficulty with the method is the fol- 
lowing. Since the quadrupole today C^(to) probes scales 
larger than the quadrupole at reionization C^(i r ) it 
seems at first that there should not be any way to make 
the method work as, strictly speaking, we are measur- 
ing Cj(i r ) which is definitely not equal to C^ito). Ex- 
amining the issue more carefully however, one sees that 
there is significant overlap between the scales spanned 
by Aj (k, to) and Aj (k,t r ) as shown in Fig.Q] Since the 
/ = 2 moment is mapped into /c-space by we expect 

to get a wide range of scales contributing to Cj(io) as 
jz(x) is broadly distributed. Moreover Aj(fc, t r ) is also 

convolved with Gj(fc) = f \/{j=U}e* T " j£° dt reTo T 
as implied by combining (|SJ| and JSJ, which further in- 
creases the overlap. The first five G;'s as well as their 
multiplication with A^(fc,t r ) are also shown in Fig. ^ 
Therefore even though what we really measure is not 
but rather the contribution to coming from the 
quadrupole at reionization, this is sufficient to reduce the 
error on significantly below the cosmic variance limit. 

One may also wonder why we could not choose some 
very early time prior to recombination to set the initial 
condition, as during tight coupling we have n(fc, t) = 
|A^(fc,f) [23j and so we can also relate the quadrupole 
with in fact the whole of the polarization spectrum, not 
just the part produced by reionization. The problem 
however in this case is that the overlap between the lo- 
cal quadrupole at the early time and A^(fc,io) would be 
minuscule. The quadrupole today will be dominated by 
the particular integral instead, and its error would be 
effectively cosmic variance again. 

To get an estimate for the error on Cj we let the am- 
plitude of A|"(fc, t r ) vary as a free parameter. This can be 
modeled by multiplying it with a parameter q by hand, 
i.e. A| , (fe, t r ) — > qA 2 (k,t r ) in the two relevant equa- 
tions and (JSJ. The variance V r of A2 (k,t r ) can then 
be estimated as V r — VqC^ where V q = Var[g] can be 
obtained using the Fisher information matrix and where 
CI = !E m (KJ 2 > = lfo°°dk fc 2 P^\A T 2 {k 7 t r )\\ As- 
suming that the only error in the polarization spectrum 
comes from cosmic variance, the Fisher matrix is a scalar 
given by 

F = E(^ + ^)i| ln ^i 2 - w 

The derivative in the expression above can be calculated 
numerically by double-sided finite difference which we 
have taken to be 5q = 0.02 around the fiducial value 
of q = 1. This is shown in FigEI The transfer functions 
were calculated using a modified version of DASh [24| . 
For the model considered we get a variance of q of 
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FIG. 2: The derivative of lnCf" with respect to q, the ampli- 
tude of the local quadrupole at the onset of reionization. 



V q = j, = 0.002. 

Since we have used the Fisher matrix to get our es- 
timate, we can assume that the posterior pdf of a r 2m 
is Gaussian, with variance V r . Therefore the posterior 
pdf of the total 02™ = eI ° Tr ^2m + ®2m wm be Gaus- 
sian with variance V2 — e^ Tr V r + Qf 7 + 2e~n> Tr C^, 
where = ±£ m <|af/J 2 > = f J dk fc 2 P*|Af'(fc)| 2 
is the variance of the particular integral and — 

^£ m R»a|£ J + c.c.) = Ijdk fc 2 p^(k,t r )A^(k) 

the variance of the correlation. Given that we have no 
extra information about the particular integral we can 
assume that its variance will not change from the cos- 
mic variance value. The variance of QT = gS m l^2m| 2 
can then be estimated as Var[C*J] = |(V2) 2 , since under 
the above assumptions, (7j will obey a x 2 distribution 
with five degrees of freedom. Our fiducial model gives 
V (QTyi^ = 0.06, which gives an error on of 24%. 

One conceptual objection with the above argument is 
that in the usual case, one has (Cj) = Cj where as 
above it is equal to V2 instead. This is not a problem as 
the true posterior pdf of C*J would no longer be \ 2 an d 
therefore its mean and variance would not be related in 
the usual way. 

As argued before, the key point is the overlap of the 
Gi{k)A^(k,t r ) with A%(k,t ). This signifies that the 
lower the reionization redshift, the better the reduction 
of cosmic variance. Since however for low reionization 
redshifts, we obtain a smaller signal in the polarization 
spectrum, we would also get a larger variance for the 
parameter q. We should therefore expect that for low 
reionization redshifts, the error on would not be im- 
proved but even be greater that cosmic variance. This 
means that as one varies the reionization redshift z r from 
zero to some large value, the error on would become 
better and better until some conspiratory value, and then 
start to become worse and worse until it reaches cosmic 
variance again. Another way to see this is the follow- 
ing. The correlation (^m^m) becomes arbitrarily small 
with increasing z r , as is implied from Fig^ which means 
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FIG. 3: Top panel : The variation of the quadrupole error 

~^tt~ (solid) and a q (dash) with reionization redshift. The er- 

ror on the quadrupole reaches a minimum as expected, around 
z = 3. The cosmic variance limit is shown as a gray line. 



Bottom panel: Variation of the ratios (solid) , ^Jr- (dash) 

C X 

and —^(dotted) with reionization redshift. 



that we are probing many independent Hubble volumes 
at that redshift hence the small variance of q. On the 



other hand 



<c 2 cj> 



is also decreasing (but slowly 



((C2) 2 ) ci 

enough) which means that the propagation of our in- 
formation on C£ to Cj becomes less effective at higher 
z r 0. This is shown in Fig.©. 

The role of other cosmological parameters is also very 
important due to imminent degeneracies, particularly 
with the optical depth. Still the quoted error above is 
a lower bound and other parameters can be included at 
a later stage along with predictions for future polariza- 
tion experiments. 

Finally let us see how a different variance on 
could affect the probabilities mentioned in the beginning. 
Strictly speaking, we need the true posterior pdf of Cj 
but based on our assumptions we can assume that it 
would be approximately Gaussian with mean given by 
the model's value and variance the value quoted above. 
Moreover, the maximum entropy principle, would also 
give a Gaussian distribution if the only knowledge about 
the distribution is the mean and the variance. If the 
fiducial model was the best fit model of a cosmic variance 
limited polarization experiment then the probability that 
the quadrupole is as low as or lower than 250 (fi K) 2 , is 
reduced to 2 x 10~ 4 . This should be taken only as an or- 
der of magnitude estimate of the actual probability which 
would have been given by the true pdf. For comparison, 
using a Gaussian distribution with variance given by cos- 
mic variance instead, we get a probability which is only 
8% different than the one quoted in the beginning. 

We have shown that it is possible to reduce the error 
on the CMB temperature quadrupole, to a value better 
than cosmic variance. The method exploits the connec- 
tion between the temperature quadrupole and the polar- 
ization spectrum generated by a period of reionization. 



This could reduce the variance of the quadrupole signif- 
icantly and has the potential to answer with confidence 
whether the quadrupole is really low or not, compared to 
a given model. The method is still at its infancy and fur- 
ther treatment is needed before it can be incorporated 
with parameter estimation techniques. It would be an 
excellent way to test new physics. 
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